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Abstract 

We have performed accurate ab-initio pseudopotential calculations for the structural and elec- 
tronic properties of ZnSe/GaAs(001) heterostructures with interface configurations accounting for 
charge neutrality prescriptions. Beside the simplest configurations with atomic interdiffusion we 
consider also some configurations characterized by As depletion and cation vacancies, motivated 
by the recent successfull growth of ZnSe/GaAs pseudomorphic structures with minimum stacking 
fault density characterized by the presence of a defected (Zn,Ga)Se alloy in the interface region. We 
find that — under particular thermodynamic conditions — some defected configurations are favoured 
with respect to undefected ones with simple anion or cation mixing, and that the calculated band 
offsets for some defected structures are compatible with those measured. Although it is not pos- 
sible to extract indications about the precise interface composition and vacancy concentration, 
our results support the experimental indication of (Zn,Ga)Se defected compounds in high-quality 
ZnSe/GaAs(001) heterojunctions with low native stacking fault density. The range of measured 
band offset suggests that different atoms at interfaces rearrange, with possible presence of vacancies, 
in such a way that not only local charges but also ionic dipoles are vanishing. 

PACS numbers: PACS: 73.40.Kp, 73.20.-r, 68.35.-p 
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I. INTRODUCTION 



II-VI/III-V heterostructures, whose prototype is ZnSe/GaAs, fabricated by molecular 
beam epitaxy (MBE) are important for electronic devices such as blue-green emitters 1 - in 
optoelectronics, spin-transistors 2 ^ and spin filters^ in spintronics. For all the proposed 
applications, it is crucial to have a control of the structural quality of these interfaces and, 
in particular, of the native stacking fault (SF) defect density which causes serious device 
degradation. 

The SFs density in II-VI/III-V interface is determined by the growth conditions. It has 
been suggested that the reduction of stacking fault density in the ZnSe/ GaAs heterostructure 
can be achieved by inserting a thin low-temperature-grown ZnSe buffer layer between the 
high-temperature-grown ZnSe epilayer and the GaAs substrate, thus obtaining SF density as 
low as ~5.4xl0 4 cm" 2 -. More recently, two procedures have been found to work in reducing 
the SF density even below 10 4 cm -2 (thus providing a very high quality system), yielding 
quantitatively similar defected densities and qualitatively similar interface compositions and 
band alignments-. 

A part from this experimental evidence, the microscopic mechanisms that control the na- 
tive defect density in II-VI/III-V structures remain still controversial. What is noteworthy 
is that from an accurate characterization of the samples with minimum stacking fault den- 
sity 6 there is experimental evidence of the formation of a ternary (Zn,Ga)Se alloy of variable 
composition with a substantial concentration of cation vacancies.— Because of vacancies, its 
average lattice parameter is smaller than the one of GaAs and ZnSe, and therefore this 
alloy is under tensile biaxial strain when epitaxially grown on GaAs substrates, accumu- 
lating a non negligible elastic energy. Evidence of formation of ordered binary defected 
compounds such as Ga2Se3 or interface layers with vacancies arranged with other symmetry 
at ZnSe/GaAs interfaces has also been reported.^ The defected compound Ga 2 Se 3 has also 
intentionally been grown epitaxially on GaAs 1 ^ and, very recently, on Si(OOl) substrates^ 
for its utilization in optoelectronics. 

The driving mechanism for the formation of this defected alloy at ZnSe/GaAs interfaces 
is still unknown and, moreover, it is not clear whether defected interfaces could be favourite 
with respect to the simplest cases of cation and/or anion mixed interfaces without vacancies. 

To this aim we have performed a comparative study based on ab-initio local-density- 
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functional pseudopotential approach of ZnSe/GaAs (001) heterostructures, including a few 
selected interface configurations with cation vacancies. Since charged interfaces are energet- 
ically unstable with respect to the interdiffusion of atoms across the interfac o 12 i 13 i 14 i 15 i 16 i 17 
we consider interfaces satisfying the charge neutrality condition. Our comparative study 
will address their structural and electronic properties as well as their formation energy. 
This work is organized as follows: in the next section we describe the theoretical and com- 
putational approach; in Sect. Ill we describe the selected interface morphologies; in Sect. 
IV and V we report our results for the structural properties and relative stability; Sect. VI 
is devoted to the discussion of the electronic properties in terms of Density of States (Sect. 
VI-A) and band alignments (Sect. VI-B); finally, in Sect. VII we draw our conclusions. 



II. THEORETICAL AND COMPUTATIONAL METHOD 

Our calculations are performed within the density functional theory framework using 
the local density approximation for the exchange-correlation functiona l 18119 with state-of- 
the-art first-principles pseudopotential self-consistent calculations. 20 The wave functions are 
expanded onto a plane-wave basis set with a kinetic energy cutoff of 20 Ry: we carefully 
check that all the relevant bulk properties of the binary ZnSe and GaAs compounds are 
well converged. The Zn-3d electrons are taken into account only via the nonlinear core 
correction which has been shown to give reliable results for ZnSe, as discussed in Ref. I"!?! 
The integration over the Brillouin zone is performed using the special k-point technique with 
a (6,6,6) Monkhorst-Pack mesh for the FCC cell and corresponding meshes for the various 
supercells. 

The theoretical lattice constant of GaAs (ZnSe) is 1.8 % (1.4 %) smaller than the ex- 
perimental lattice constant a^, As =5.65 A3i (af x " 5e =5.67 A 21 ) but the small relative lattice 
mismatch is well reproduced (< 0.8 %). In our calculations, we neglect this tiny lattice 
mismatch and we fixed the in-plane lattice constant to the theoretical one of GaAs which 
is typically the substrate in experimental samples. The ionic degrees of freedom are fully 
taken into account by optimizing the atomic positions via a total-energy and atomic-force 
minimization—, with a threeshold of 1 mRy a.u. -1 for atomic forces and of 0.01 mRy for 
total energy changes between two consecutive relaxation steps. 

The interfaces are modeled by tetragonal supercells with periodic boundary conditions 
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and contain a slab for each constituent material. We study the relative stability by consid- 
ering the interface formation energy:— 

2£jintf {-^supercell L~ii=l "r i q\ 

' ran 

where E superceU is the calculated total energy of the supercell, N species is the number of the 
chemical species involved (which are 4 in the present case), ri 1 the number of atoms of the 
species i and fi l is the corresponding chemical potential; (m x n) is the reconstruction (the 
in-plane periodicity of the interface), so that N = mn is the number of atoms in the two- 
dimensional interface unit cell and Ej 1 *^ refers to a (lxl) interface area. In general E™ 1 ^ 
refers to a mean value of the formation energy of the two possibly inequivalent interfaces 
present in each supercell. We calculated the bulk chemical potentials considering the ele- 
mental forms of Zn (hep),— Se (trigonal),— Ga (orthorhombic),— As (trigonal),— GaAs and 
ZnSe (cubic). 21 We employed the Gaussian smearing technique 2 ^ for Brillouin integration 
of metallic systems. For each bulk elemental form we have chosen the proper fc-point mesh 
which gives well converged values of the corresponding chemical potential, with a numerical 
uncertainty less than ^20 meV per bulk formula units. The calculated heat of formation 
of bulk ZnSe and GaAs are —1.61 eV and —0.80 eV respectively which well compare with 
the experimental values (—1.67 eV and —0.84 eV) ! 21 ^ this indirectly tests the reliability 
of our calculated elemental chemical potentials. They also agree with previous theoretical 
calculations . MM 



III. INTERFACE STRUCTURES 



We discuss in this Section the selected interface morphologies that we use as simplest 
models to describe the defected ZnSe/GaAs (001) junctions. Based on experimental sug- 
gestions about the formation of (Zn,Ga)Se compound in a As-depleted region^ and for the 
sake of simplicity we consider only cation vacancies and we exclude the possibility of other 
defects such as interstitials or antisites. 

Simple electrostatic considerations, already pointed out by Harrison along time ago,^ 
suggest that no ionic charge accumulation can occur at the junction since uncompensated 
charges would set up an electric field extending throughout the overlayer, with an energy 
accumulation clearly unfavourable as long as the overlayer thickness increases.— The charge 
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neutrality condition at the interface could be expressed in terms of ionic charges by cal- 
culating their macroscopic average density profile along the direction of the junction and 
imposing that possible deviations are compensated. 32 Alternatively, it can be formulated 
also in terms of bond charge compensation.— In general, at interfaces between heterovalent 
constituents there are "wrong" chemical bonds with either more ("donor bond") or less 
( "acceptor bond" ) than two electrons.— A practical rule for the compensation of donor and 
acceptor bonds can be derived as follows. In the zincblende structure between two consecu- 
tive atomic planes % and % + 1 in the [001] direction there are two bonds for an area of cLq/2, 
where ao is the bulk lattice constant, with two electrons each one and hence with a total of 
four electrons. These electrons come from the contributions Zi/2 and Z i+1 /2, where Z ; is the 
ionic charge of the atoms in the i plane. Zi and Z i+1 are 3 and 5(2 and 6) in GaAs(ZnSe) 
and therefore trivially 

(Z M+i ) = (Z i + Z i+1 )/2 = 4, (2) 

where (Z iji+1 ) is the average of bond electrons between the two corresponding consecutive 
atomic planes. The compensation of donor and acceptor bonds requires that Eq. |2] holds 
also in the interplanar region at the interface, with Z; referring in general to the average 
ionic charge per atom in the i plane. If not satisfied over a pair of consecutive atomic planes, 
the condition must be fullfilled by compensation with neighboring bonds, i.e. in neighboring 
pairs of atomic planes. 

We first consider, for the sake of example and comparison, the simplest interface mor- 
phologies with atomic intermixing limited to one or two atomic planes, without vacancies 
and satisfying the charge neutrality, already addressed in Refs. ^ 12- The composition 
profiles for these simple cases are: 

1C: 1-plane Cation-mixed: • • • -Ga-As-(Gai Zni)-Se-Zn-- • • 

v 2 2 ' 

1A: 1-plane Anion-mixed: ■ ■ • -As-Ga-(Asi Sei)-Zn-Se-- ■ • 

v 2 2 ' 

2CA-GaSe: 2-planes Cation- Anion mixed, Ga and Se rich: 

• • • -Ga-As-(Zni Ga3)-(Se3 Asi)-Zn-Se-- • • 

4 4 4 4 

2CA-ZnAs: 2-planes Cation-Anion mixed, Zn and As rich: 

• • -As-Ga-fSei As3)-(Zn3 Gai)-Se-Zn-- • • 

v 4 4 ' v 4 4 ' 

The profiles of the average ionic charges Z; on the atomic planes are: 



1C: 3-5-2.5-6-2 

1A: 5-3-5.5-2-6 

2CA-GaSe: 3 - 5 - 2.75 - 5.75 - 2 - 6 

2CA-ZnAs: 5 - 3 - 5.25 - 2.25 - 6 - 2 

and the sequences of the average bond electrons in the interplanar spaces (Z ii+1 ) are: 

1C: 4- 3.75 -4.25-4 

1A: 4-4.25 - 3.75-4 

2CA-GaSe: 4 - 3.875 - 4.25 - 3.875 - 4 

2CA-ZnAs: 4-4.125- 3.75 -4.125 -4 

It is easy to check that for all these interfaces Eq. |21 is satisfied by a compensation 
occurring at most over three interplanar spacings. The latter interfaces have no ionic dipole, 
as it would be the case for the abrupt, nonpolar (110) interface. The interfaces 1C, 1A and 
2CA-ZnAs, 2CA-GaSe are related by a simultaneous exchange of Ga with As and Zn with 
Se. 

Introducing cation vacancies, many other different morphologies compatible with the 
charge neutrality condition are possible. We generalize the case of 1-plane Cation-mixed 
interface by considering the possibility of another Cation-mixed plane with vacancies as 
follows: 

• • • — Ga — As — (Ga x Zn y Vi- x ^y) — Se — {Ga w Zn z V\- w - z ) — Se — Zn — ■ ■ ■ 

where V indicate vacancies and x, y, w, z are the in-plane atomic concentrations of the 
corresponding atomic species, with the obvious condition: 

x + y < 1 and w + z < 1. (3) 

For this case, the sequence of the average ionic charges Zi on the atomic planes is: 

3 - 5 - (3a: + 2y) - 6 - (3w + 2z) - 6 - 2 



6 



and the corresponding sequence of (Z i i+ i) as defined above is: 

| 5 + 3a; + 2?/ 3x + 2y + 6 6 + 3w + 2z 3w + 2z + 6 j 
2 2 2 2 

The charge neutrality condition, that is Eq. applied to this composition profile gives: 

Q(x + w) + A(y + z) = 9 (4) 

We note that Eq. H]is symmetric under the exchange of x <-> w and/or y <-> z. In the 
following we focus only on the pairs of complementary solutions related by both x <-> w and 
y «-> 2 exchange and satisfy Eq. El this corresponds to a swap of the cation planes across 
the Se plane. 

First, we are going to consider systems with a high (H) local concentration of vacancies: 
in particular an interface plane with 50% concentration of vacancies, a configuration not 
very realistic but that can be easily simulated with small supercells. This case corresponds 
tox=^,w = l,y = z = and the complementary solution, with the following composition 
profiles: 

1HV: ■ ■ ■ — Ga — As — Gai — Se — Ga — Se — Zn — ■ ■ ■ 

2 

lHV-swap: ■ ■ ■ — Ga — As — Ga — Se — Gai — Se — Zn — ■ ■ ■ 

2 

These profiles are schematically shown in Fig. 1 together with a sketch of the bare 
ionic charge profile, averaged over two consecutive atomic planes and normalized to a bulk 
zincblende cell, indicating the extension and the value of the interface dipole. 

Next, we consider the solutions with z—1 and w=0 which reduce to the case of vacancies 
just confined in a single plane between GaAs and ZnSe. Out of the possible cases, we select 
the one with x = | and y = 0. Considering also the swap across the Se plane, we have the 
following composition profiles: 

IV: ■ ■ ■ — Ga — As — Ga^ — Se — Zn — Se — Zn — ■ ■ ■ 

6 

1 V-swap: ■ ■ ■ — Ga — As — Zn — Se — Gas — Se — Zn — ■ ■ ■ 

6 

In terms of ionic charges, each Zn is equivalent to | of Ga, so that if we substitute Zn with 
Ga2 in the composition profiles of the possible neutral interfaces we still maintain the charge 
neutrality condition and we have: 
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2V: ■ ■ ■ — Ga — As — Gas — Se — Ga-i — Se — Zn — ■ ■ • obtained from IV; 

6 3 

2V-swap: • • • — Ga — As — Ga2 — Se — Gas — Se — Zn — ■ ■ ■ obtained from lV-swap. 

These interfaces have two planes with vacancies and correspond tox = ^, w = ^, y = z = 
and the complementary case x = ^, w = ^, y = z = 0. 

All the interfaces considered above are characterized by the presence of Ga-Se and/or 
Zn-As bonds, which are not present in the bulk constituents. Considering that recent ex- 
perimental results do not indicate the presence of Zn-As bonds^ we focus here mainly on 
structures with Ga-Se bonds and we include for comparison only one with Zn-As bonds, 
namely the one labelled as "lV-swap". Remarkably, in the entire range of the chemical 
potentials involved, the interface having Zn-As bonds is thermodynamically unstable with 
respect to the complementary one (IV) which has the same stoichiometry but no Zn-As 
bonds (see Sector©. 

IV. STRUCTURAL PROPERTIES 

Because of the cation vacancies, lattice distortions are sizeable and affect mainly the anion 
sublattice. Below, we briefly discuss the structural properties of the different configurations 
studied. Lattice distortions have non negligible effects on electronic properties and stability 
of the system, which will be discussed in the next Section. 

A. Interfaces with 50% vacancies layer 

Both 1HV and lHV-swap supercells are characterized by a huge local concentration of 
cation vacancies, having a plane with 50% of vacancies. The supercells considered have 2 
atoms per layer and 20 atomic layers (the total length of the supercell is 33.27 A). Atomic 
relaxations are important in the interface region, and result in sizeable variations of the bond 
lengths — involving in particular of acceptor and donor bonds — and interplanar distances 
along the growth direction, which are reported in Fig. 2. 
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B. Interfaces with different vacancy concentration 

The supercell describing the 1(2)V and l(2)V-swap interfaces have 6 atoms per layer 
and 20 atomic layers, with a total length of the supercells of 27.72 A. In Fig. [2] we show 
the average interplanar distances along the [001] direction. In this cases, we have a non 
negligible buckling of the atomic planes, in particular in the interface region, which reduces 
progressively far from the interface, although propagating also several atomic layers far away, 
much more in 2V and 2V-swap than in IV and lV-swap. In Fig. 2 we report the average 
interplanar distances with an error bar determined by the maximum and minimum interlayer 
distances between adjacent atomic layers. 

Although the composition profile along the growth direction is symmetric with respect 
to a middle plane in the constituent slabs, the two interfaces in the simulation cell are not 
equivalent in their 3D structure, and this justifies the lack of symmetry in the pattern of the 
interatomic distances. Incidentally, we note that the larger is the buckling in the interface 
region, the stronger is the effect propagating in the bulk slabs. 

V. THERMODYNAMIC STABILITY OF INTERFACES 

We study the relative stability of the interfaces using Eq. (JTJ). Although the precise values 
of the chemical potentials are unknown, strongly depending on the growth process and on 
the local environment, their range of variation and mutual relationships can be established 
in condition of thermodynamic equilibrium. In particular: 

VgIas = ^Ga + ^As and [l b ZnSe = ^Zn + VSe (5) 

for the equilibrium between the interface and the bulks; 

IH < f4 ulk (6) 

to exclude the formation of precipitates (since when ^ = ^ ulk the gas phase condensates to 
form the elemental bulk phase). We could also exclude the formation of other precipitates, 
like Zn3As2 and Ga2Se3, but this would only further restrict the possible range of variations 
of the chemical potentials involved, without changing the main conclusions of our study. 
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By defining the heat of formation of AB compound as: AHab = I^ab ~ ^A k ~ $B k ' ' 



so that a negative AHab means exotermic reaction, one gets the following relations: 

AH GaAs + /^f < Has < (7) 

and 

AH ZnSe + ^ < < pg* (8) 
We choose \iz n — ^Ga and fise as basic variables, and with some algebra we obtain: 

/4f - »Ga + ^H ZnSe < VZn ~ VGa < ~ »Ga ~ ^H GaAs (9) 

We can introduce new variables 

fi Zn ~ JlGa = liZn ~ ^Ga - (^n ~ (10) 

and 

fa = »Se - 1%!* (11) 

and discuss the relative stability of the interfaces in the ranges of variations of these variables. 
We can specify Eq. (JTJ) as follows: 



2J5«**/ _ Esupercell ~ Nse^ZnSe ~ ^As^GaAs + {^Se ~ N Zn )H>Zn + {N As ~ N Ga )fI Ga 



The supercells with no vacancies that we have considered here are stoichiometric 
(N se=N zn=N znSe, and NGa=NAs=NG aJ 4s), therefore their formation energy is independent 
on the chemical potentials and Eq. (fT2"|) reduces to a constant value. The resulting interface 
formation energies for 1C and 1A are nearly degenerate, within 10 meV per (lxl) interface 
unit cell. The same holds for 2CA-ZnAs and 2CA-GaSe. The relaxations lower the forma- 
tion energies by 60 meV (50 meV) per interface unit cell for 1C and 1A, whereas the effect 
is stronger for 2CA-ZnAs and 2CA-GaSe where the formation energy is lowered by ~ 120 
meV per (lxl) interface unit cell. 
For the other structures, from Eq. (fT2*|) we obtain: 

E f!lHV = C 1HV + (P<Zn - J^Ga) ~ (13) 



E jlv = C IV + \{^Zn ~ pGa) ~ ^Se (14) 



2' 
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Ef*2V ~ ^2V + (HZn — UGa) — (15) 

We observe that structures related by a swap of atomic planes have the same iVj's, so they 
have the same expressions for the interface formation energy: therefore Eqs. (|1H|14[ ITK|) also 
hold for the complementary structures. The constants are characteristic of the specific 
interface and do not depend on the length of the supercell used: test calculations show that 
they change only by ~ 10 meV per (lxl) interface unit cell when increasing the length of 
the supercells: this is a further indication that the size of our supercells is large enough. 

As it can be seen from Eqs. (|1H[ I14|15j) . the interface energy Ej 1 '^ is a linear function 
of fizn — fi-Ga and Jlse with coefficients determined by the particular stoichiometry. It is 
convenient to discuss the dependence of Ej n ^ on an individual variable separately, setting 
the other one equal to its limiting value given by Eqs. (|8|9[) . In Fig. El we show the interface 
formation energies as a function of p, Zn ~ fi>Ga for lower (upper) limit of p,s e (top) and as a 
function of Jxs e for lower (upper) limit of Jl Zn ~ fica (bottom). The horizontal line corre- 
sponds to interfaces with no vacancies, namely 1C, 1A, 2CA-ZnAs, 2CA-GaSe which appear 
degenerate on the scale used. Right (left) hatched areas represent the range of variation 
of formation energy for 1HV, lHV-swap, 2V, 2V-swap (IV, lV-swap) configurations. The 
lowest formation energy for right (left) hatched area corresponds to 1HV (IV) configuration. 

Now we discuss the limiting cases corresponding to Jige = AH Zn se and JLs e = 0, i.e. the 
lowest and uppermost values of fls e (see Fig. El top part). The heat of formation for ZnSe 
determines the lower limit of Jis e (see Eqs. I8|lljl . Ej 1 '^ is a function of the other chemical 
potential. It can be seen that the most stable structures remain the same (1HV and the 
mixed ones with no vacancies) independently on the value of Jls e , although the value of 
(Uzn ~ Ugo) at which they have competing energies changes a little bit. As a general trend, 
high values of p,s e should stabilize the formation of defected interfaces (in Fig. El top part, 
compare the M&x{jlse} respect to Min{/is e } case). This observation agrees with previous 
calculation dealing with GaAs(OOl) surfaces where it has been shown that for an increasing 
Se chemical potential (high limit of Jlse) the formation of a Ga vacancy beneath the surface 
becomes energetically favourable, driving the surface stoichiometry towards Ga2Se3.— On 
the other hand, substantial changes occur considering the extreme cases for Ji Zn — J^Ga (see 
Fig. El bottom part). This indicates that the variation of Jlse has much less effect than a 
variation of jl Zn — JiGa- This can be also argued from the different relative weight of the 
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basic variables in Eqs. [T31 EH UB 

In the high Jige limit (see Fig. El top part), the defected interfaces are favored over 
the mixed ones for almost the whole range of variation of p,z n ~ J^Ga- At variance, in the 
limit of high jlz n ~ fica (see Fig. El bottom part) all the defected interfaces turn out to be 
unstable over the undetected mixed-ones, irrespectively on the variation of p,g e ] in the limit 
of low p,zn ~ fiGa, most defected interfaces are favored over the mixed-ones. Therefore a low 
fizn — fi>Ga growth condition should favour the formation of defected interfaces. 

For comparison with experiments, we observe that ZnSe is usually grown over a GaAs 
buffer, in a MBE chamber. Therefore the physically independent variables are fizn and 
Use- As far as fige is concerned, the uppermost value corresponds to Se-rich condition: 
in case of evaporation in UHV, like in the MBE chamber, this would be achieved with a 
low Zn/Se beam pressure ratio and consequently a high Se flux, with possible formation 
of cluster with a bulk-like crystal structure from the Se atoms deposited on surface. The 
physical interpretation of the other variable, fizn~ ^Ga, is less straightforward. Moreover, in 
a MBE chamber, \iQ a can be considered as a constant as only Zn and Se atoms come from 
the effusion cells. Therefore the uppermost value of /i^ n — \iQ a should correspond to rich 
Zn condition (and viceversa). This interpretation in consistent with the above discussion 
about the dependence of formation energy on the basic variables. We also point out that, 
in principle, a different choice of basic variables is possible but this would not have affected 
the main results. Our choice is the simplest one. 

We observe that 2V, 2V-swap, 1HV, lHV-swap configurations have the same linear de- 
pendence on Jlzn — V-Ga and Jlge so that we can easily compare their relative stability in the 
whole range of variability of thermodynamic variables. Their relative order, from the most 
stable, is: 1HV, 2V, 2V-swap, lHV-swap. If we define an average number of "wrong bonds" 
(including the fictitious bonds defined with the vacancies) per anion in the supercell, we can 
empirically observe that (at least for the relative stability among these structures) the most 
stable interface is the one with the lowest number of Se-Ga and As-V bonds and the highest 
number of Se-V bonds. On the other hand, IV and lV-swap are stoichiometric but differ for 
the presence of Zn-As bonds in the latter. As mentioned before, the lV-swap configuration 
is unstable with respect to IV, in the whole range of chemical potentials. This result is 
consistent with the experimental indication of Ga-Se bonds but no Zn-As bonds and with 



the results of Ref. 



35 which is however limited to undetected interfaces and does not take 
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into account the charge neutrality prescription. 

Finally we caution the reader that we have just considered trends in relative stability 
among different interfaces only on the basis of thermodynamic arguments, and therefore 
the high relative stability of the interface with unrealistic concentration of vacancies has 
not an absolute validity. Kinetic effects should play an important role on stability and in 
particular could hinder the formation of interfaces with unrealistic vacancy concentration 
(see also discussion in Section IVIB|) . Taking into account kinetic effects would go beyond 
the scope of the present study. 

VI. ELECTRONIC PROPERTIES 
A. Density of States 

The presence of vacancies induces in general strong perturbations on the electronic states 
due to the large fluctuations of the crystalline field and may induce localized states at 
different energies. As an example, we report in Fig. 4 the total Density Of States (DOS) 
and the atomic Projected Density Of States (PDOS) on different atomic layers for the 1HV 
and lHV-swap structures. In particular, the PDOS is shown from the topmost to the bottom 
panel for the sequence of atomic planes along the growth direction from the ZnSe side to the 
GaAs side through the plane with vacancies. It is also shown the "excess" DOS (black area), 
defined as the difference between the PDOS and the corresponding bulk DOS, if positive. 
The vertical arrows in each panel denote the energy position of the cation and anion s states. 

From the total DOS (topmost panels) it is evident that the main difference between the 
two structures with a high vacancy concentration in the interface region is the density of 
states in the forbidden energy gap: the structure 1HV is metallic, whereas the lHV-swap 
structure maintains the semiconducting character of its constituents. Atomic relaxations 
tend to reduce the effects and to clean the gap, but not completely, as it can be seen 
comparing the two topmost panels for both 1HV and lHV-swap structures (solid/dashed 
line for relaxed/unrelaxed structure) 

Looking at the PDOS, it is clear that the gap states in the total DOS of the 1HV structure 
originate mainly from states localized in the interface region where the atomic species have 
"wrong" chemical bonds, with the major contribution deriving from the Ga layer sandwiched 
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between two Se layers and to the Se atoms bonded to Ga atoms and vacancies. An effect 
of the strong local fluctuations of the crystalline potential in the interface regions is the 
variation in the energy position of the s-Se peaks with respect to their value in the bulk 
region and variation in the shape of the PDOS of Ga atoms bond to Se and As (see the 
arrows) . 

By inspection of the PDOS in the bulk sides we can see a large difference in the relative 
positions of the bulk s-Se and s-As peaks across the interface in the 1HV and lHV-swap 
structures: in 1HV, the s-Se peaks are lower with respect to the s-As peaks by about 
4 eV, whereas in the lHV-swap structure this energy difference is strongly reduced and 
almost vanishing. This is related to the different ionic interface dipole due to the atoms 
swap, sketched in Fig. 1, and results also in different bands alignments, as discussed in the 
following Subsection. 

Some of the effects of the fluctuations of the crystalline field can be explained with a 
simple model accounting for the "wrong" chemical bonds. We characterize these "wrong" 
bonds by defining the excess average ionic charge of the nearest-neighbors (NN), Aznn^ 
averaging also over the vacancies: if Az^n is positive (negative), the higher (lower) it is, 
the more (less) attractive it is the local potential with respect to a bulk environment. 

In 1HV structure, the As atoms close to the plane of vacancies have Aznn=— 0.75: they 
are in a crystalline potential less attractive than in the bulk and their valence band states 
are shifted upwards in the direction of the forbidden energy gap. The other atoms in the 
interface region have a positive Az^n' Aznn = 0.5, 1.0, 0.25, 0.5, respectively, for the Se 
layer between Zn and Ga, for Ga between Se layers, for Se between Ga layers, and for the 
Ga plane with vacancies. A crystalline potential more attractive than the bulk one pushes 
down in energy the states from the conduction band. This roughly explains the origin of 
the gap states in this structure. 

In the lHV-swap structure, the local fluctuations of the crystalline field are weaker. 
Se atoms between the Zn and the Gao.5 layer with vacancies show electronic states not 
very different with respect to the bulk; and Zn atoms nearest to the interface region (with 
Az/vat=0) also show features similar to the bulk. The Se atoms sandwiched between two Ga 
planes (with Az^^=0.25), strongly differ from the bulk case but not in the region of the 
energy gap. 

Extending our analysis to the other defected interfaces studied here, we notice that the 



14 



presence of electronic states within the energy gap is a common feature, and the lHV-swap 
structure is an exception. 



B. Band alignments 

We summarize schematically in Fig. |3] the band alignments calculated for the selected 

321 . splitting the band offset 



structures considered here. We follow the approach of Ref. 
into two contributions: V BO = AE V + AV, where AE V is the band structure term, i.e. the 
energy difference between the relevant valence band top edges of the two materials measured 
with respect to the average electrostatic potential in the corresponding bulk crystal and AV 
is the electrostatic potential lineup containing all interface-specific effects and extracted from 
supercell calculations. The potential lineup in turn includes a contribution related to the 
ionic dipole and its screening (indicated as AV hetero in Ref. I32T ) plus a term which is purely 
electronic and independent on interface details. The spin orbit effects are added a posteriori 
using experimental data. 

As predicted for the heterovalent heterostructures, the band alignments are strongly 
dependent on the specific interface morphology. This is particularly true in the occurrence of 
vacancies, as it can be easily seen by inspection of Fig. 5: for undefected (without vacancies) 
interfaces the maximum variation is of the order of 1 eV (compare 1C and 1A cases) , 15 i 37 i 38 
whereas it is much larger in case of vacancies by changing their concentration and position, 
the extreme variation (of about 4 eV) being found for 1HV and lHV-swap structures, where 
the alignment changes from the so-called "broken gap" type (1HV) to "staggered" (1HV- 
swap). We point out that such a huge variation is associated to unrealistically high local 
concentration of vacancies. In real samples we expect vacancies diluted over a larger interface 
region. 

In the spirit of linear response theory (LRT) 12 i 32 ' 3 ? the variation of band offset in heterova- 
lent heteroj unctions has been rationalized and quantitatively explained in terms of the line- 
ups associated to the different interface ionic charge distribution and the corresponding elec- 
tronic screening, i.e. of AVhetero- Neglecting details such as optimized atomic positions, sim- 
ple calculations based on elementary electrostatics and LRT give AV/ ie t ero =7re 2 /2ao(e) ~ 0.5 
eV for 1C, IV and 2V, the opposite for 1A, AV/ ietero =0 for 2CA-ZnSe and 2CA-GaSe, and 
AV hetero =9ne 2 /2a (e) « 4 eV for 1HV and AV /ietero =-3vre 2 /2a (e) « -0.33 eV for 1HV- 
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swap, where ao is the lattice constant and (e) a proper average^ between the dielectric 
constant of GaAs and ZnSe (which are 10.86 and 5.7 for GaAs and ZnSe respectively^). 
By comparing the LRT predictions with the full self-consistent calculations reported in Fig. 
5 we can comment about the validity of the LRT: its predictions are quantitatively correct 
for undefected structures, but fails with discrepancies up to 1 eV in case of vacancies (LRT 
predicts band offset IV and 2V equal to 1C, and lV-swap equal to 2V-swap). It is not 
surprising: the LRT better holds in case of interface configurations without vacancies since 
the perturbation due to the presence of the interface is related to rather small differences 
between atoms (e.g. Zn vs. Ga, Se vs. As), and it is expected to be less accurate in case of 
strong perturbations such as vacancies. 

The measured variation of the band offset in real samples obtained with different growth 
procedures is much lower than what found here for different vacancy configurations: it is 
limited to about 0.6 eV, as indicated by the range of the experimental data in the right hand 
side of Fig. 5 with the label "Exp(2)". We show separately the right hand side of Fig. 5 
with the label "Exp(l)" the experimental values from Ref. [D, ranging from 0.60±0.05 eV to 
0.72±0.10 eV for the two samples with lowest SF density, whose composition of the interface 
region is indicated to be 45 at% Se, 34 at% Zn, 21 at% Ga and 50 at% Se, 40 at% Zn, 10 
at% Ga respectively. A common indication extracted from the experiments is a type-I or 
"straddling" band alignment (the band edges in the GaAs slab are both within the band 
gap of ZnSe): remarkably, this is qualitatively compatible with most configurations studied 
here, including some with vacancies. 

It is not possible to extract any precise indication on the morphology of real samples from 
the comparison between experimental results and numerical predictions: there are many 
different configurations which can give the same band offsets and, moreover, relative stability 
crucially depends of the kinetics of the growth process. However, we can conclude that the 
occurrence of vacancies is compatible with the experimentally measured band offsets, and 
the most likely configurations in real samples seems to be characterized not only by charge 
neutrality but also by vanishing ionic dipoles. 
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VII. CONCLUSION 



We have studied with accurate ab-initio calculations several selected interface configu- 
rations for the ZnSe/GaAs(001) heteroj unctions satisfying the charge neutrality condition, 
including also defected configurations with cation vacancies. We have shown that in some 
particular thermodynamic conditions (basically rich-Se chemical potential) the formation 
of defected interfaces with cation vacancies is even favoured over the simpler, undetected, 
unstrained, anion- and/or cation-mixed interfaces. We would like to note that one of the 
procedure leading to interface with low SF's density (and vacancy at interface) involve fab- 
rication of ZnSe buffer layer at interface in highly Se-rich condition.— The band alignments 
are strongly dependent on the particular interface morphology, in particular in the occur- 
rence of vacancies; on the other hand, different morphologies can correspond to the same or 
to similar band offsets, so that a unique correspondence between morphology, band offset, 
relative thermodynamic stability cannot be established. Nevertheless, the predicted sta- 
bility of some defected interfaces and the compatibility of predicted band alignments with 
measurements support the experimental evidence of (Zn,Ga)Se defected compounds in high- 
quality ZnSe/GaAs(001) heteroj unctions with low native stacking fault density. The range 
of measured band offset suggests that different atoms at interfaces rearrange, with possible 
presence of vacancies, in such a way that not only local charges but also ionic dipoles are 
vanishing. 
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FIG. 1: Composition and bare ionic charge profiles (in unit of electronic charge per zincblende 
unit cell) in the interface region of the defected heterostructures 1HV, lHV-swap, IV, lV-swap, 
2V, 2V-swap along the (001) growth direction. The ionic charge is averaged over pairs of adjacent 
(001) atomic layers and normalized to the bulk zincblende unit cell. Interface cationic planes are 
emphasized in boldface. 

FIG. 2: Interplanar distances (in A) for the different structures considered including cation vacan- 
cies. The average value of the interplanar distances are indicated by closed symbols, and the error 
bar indicate the range of local variation due to the possible buckling of the atomic planes. 

FIG. 3: Top: Formation energy [in eV/lxl surface] of the mixed and defected interfaces as a 
function of the difference of the Zn and Ga chemical potentials, for lower(upper) limit of p, Se . Dotted 
lines refer to the undetected 1C, 1A, 2CA-ZnAs, 2CA-GaSe interfaces, which appear degenerate on 
this energy scale. Hatched areas indicate the range of interface formation energy for interfaces with 
the same stoichiometry: right (left) hatched areas represent the range of variation of formation 
energy for 1HV, lHV-swap, 2V, 2V-swap (IV, lV-swap). Those with the lowest formation energy 
are lHV-swap and IV for the two families respectively. Bottom: as before, as a function of the 
difference of the Se chemical potential for the lower (upper) limit of Ji Zn — Jl Ga . The chemical 
potentials have been rescaled with respect to their correspondent bulk value (see text). 

FIG. 4: Total DOS (upper panels) and atomic Projected Density of States over layers along the 
growth direction for the 1HV (panels a) and lHV-swap (panels b) interface. The "excess" DOS is 
shown as black area (see text). Arrows indicate the energy position of the cation and anion s states 
in each layer. Energies are referred to the Fermi energy for 1HV and to the Valence Band Maximum 
for lHV-swap. Solid/dashed line in the topmost panels are for relaxed/unrelaxed structure. 
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FIG. 5: Scheme of the band alignments for the different interfaces investigated. The relative energy 
scale is referred to the band edges of GaAs, which are indicated on the left, and the relative position 
of ZnSe is then indicated for the different morphologies. The horizontal dotted lines refer to GaAs 
band edges. The label "Exp(l)" refers to experimental values in Ref|3 and "Exp(2)" to all the 
other measurements reported in the literature-^. 
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